Calculation of Magnetic Field Noise from High-Permeability Magnetic Shields and 

Conducting Objects with Simple Geometry 
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High-permeability magnetic shields generate magnetic field noise that can limit the sensitivity 
of modern precision measurements. We show that calculations based on the fluctuation-dissipation 
theorem allow quantitative evaluation of magnetic field noise, either from current or magnetization 
fluctuations, inside enclosures made of high-permeability materials. Explicit analytical formulas 
for the noise are derived for a few axially symmetric geometries, which are compared with results 
of numerical finite element analysis. Comparison is made between noises caused by current and 
magnetization fluctuations inside a high-permeability shield and also between current-fluctuation- 
induced noises inside magnetic and non-magnetic conducting shells. A simple model is suggested 
to predict power-law decay of noise spectra beyond quasi-static regime. Our results can be used to 
assess noise from existing shields and to guide design of new shields for precision measurements. 



I. INTRODUCTION 

Passive magnetic shields are frequently used in preci- 
sion measurements to create a region in space which is 
magnetically isolated from the surroundings [l[. A few 
layers of nested shells made of high-permeability metals, 
such as mu-metal, routinely provide in table-top exper- 
iments a quasi-static shielding factor in excess of 10 4 . 
Such a shield, on the other hand, generates thermal mag- 
netic field noise which often exceeds the intrinsic noise of 
modern high-sensitivity detectors such as superconduct- 
ing quantum interference devices (SQUIDs) and high- 
density alkali atomic magnetometers [2j. 

Magnetic field noise generated by thermal motion of 
electrons (Johnson noise current) in metals has been 
much studied in the past in the context of applications 
of SQUID magnetometers [U, Q, and more recently as a 
source of decoherence in atoms trapped near a metallic 
surface Q. A majority of these works were devoted to 
low frequency noise from Johnson noise current in non- 
magnetic metals. A few authors also considered noise 
from high-permeability metals of flat geometry The cal- 
culations presented in these works, however, were not 
particularly amenable to extension to other geometries, 
such as those of cylindrical shields often used in table-top 
experiments. Nenonen et al., for example, used calcula- 
tion of noise from an infinite slab to estimate noise in- 
side a cubic magnetically shielded room for biomagnetic 
measurements [6|. As shown below, the validity of such 
extrapolation is not immediately clear, given the image 
effect of high-permeability plates. Lack of explicit for- 
mulas and qualitative scaling relations for magnetic field 
noise from high-permeability shields have caused some 
confusion about the contribution of such noise in certain 
experiments (See discussions in Refs. 

Among different strategies that have been demon- 
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strated to calculate magnetic field noise [1,0, Q, a partic- 
ularly versatile method is the one based on t he g eneral- 
ized Nyquist relation by Callen and Welton 0, which 
later led to the fluctuation-dissipation theorem. Here the 
noise from a dissipative material is obtained from calcu- 
lation of power loss incurred in the material by a driving 
magnetic field. Sidles et al., for example, presented a 
comprehensive analysis of the spectrum of magnetic field 
noise from magnetic and non-magnetic infinite slabs with 
a finite thickness using this principle [ill ]. 

A particularly useful feature of the power-loss based 
noise calculation is that it allows calculation of noise 
from multiple physical origins, including Johnson noise 
current in metals and domain fluctuations in magnetic 
materials. The noise of the latter kind in ferromag- 
nets, which can be associated with magnetic hysteresis 
loss, was previously studied for toroidal transformer cores 
where field lines were confined in the core material [T3 |. 
In a recent work [l3[ Kornack et al. measured magnetic 
field noise in the interior of a ferrite enclosure with an 
atomic magnetometer which was consistent with predic- 
tions based on numerical calculation of power loss in the 
ferrite. The same paper also presented results of ana- 
lytical calculations of the noise inside an infinitely long, 
high-permeability cylindrical tube. 

In this work we show how similar calculations can be 
performed for other geometries with cylindrical symme- 
try, and derive a general relationship between magnetic 
field noises from current and magnetization fluctuations 
in shields with such geometries. For metallic shields, we 
show that the Johnson current-induced noise is cither 
suppressed or amplified, depending on the shape of the 
shield, due to a high permeability. This partly explains 
previous confusion about noise contributed by magnetic 
metals. Analytical calculations leading to our key results 
were confirmed by numerical calculations on representa- 
tive geometries using commercial finite element analysis 
software. In order to explain frequency dependence of 
noise from metallic and magnetic plates reported in lit- 
erature, we propose a simple model which correctly pre- 
dicts observed power-law decays in noise spectra. We also 
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present in the Appendix analytical calculations of noise 
from non-magnetic conducting objects that can model 
other common experimental parts used in precision mea- 
surements. 



II. PRINCIPLES 

The principle of calculating magnetic field noise from 
energy dissipation in the source material has been demon- 
strated by several authors. For example see Refs. [!, H, 
HTj . The arguement is summarized as follows. If at a 
point f there is a fluctuation of magnetic field along di- 
rection n, given by its power spectral density Sb (/), an 
N-turn pickup coil located at r directed along n will de- 
velop a fluctuating voltage, according to the Faraday's 
law, with power spectral density 



S v (/) = AWuPSb (/) 



(1) 



Here u> = 2irf and A is the area of the pickup coil, as- 
sumed to be small so that the field is uniform over the 
area. We further assume that the coil is purely induc- 
tive, for example by making it superconducting, so that 
in the absence of an external material (noise source) there 
is no voltage fluctuation due to conventional Nyquist 
noise, Sv,coU — ^kTR co u = 0. Now assume that we 
take the pickup coil and the material responsible for 
the noise as a single effective electronic element, whose 
small-excitation response is characterized by an impe- 
dence Z. The fluctuation-dissipation theorem applied to 
this system states that the voltage fluctuation at the ter- 
minals of the pickup coil is related to the real part of Z, 
Re[Z(f)} = R eSl by 



S v (f)=AkTR eS (f) 



(2) 



Here the system is assumed to be at thermal equilib- 
rium at temperature T, k is the Boltzmann constant, and 
the effective resistance P c ff is obtained from the time- 
averaged power dissipation in the system 



P{f)= l -I 2 R cS (j) 



(3) 



incurred by an oscillating current I (t) = I sinwt flowing 
in the pickup coil whose amplitude / is small so that the 
response is linear. In the absence of the resistance of the 
pickup coil itself, the power dissipation is entirely due to 
the loss in the material driven electromagnetically by the 
current / (t) . From Eqs. (|1I2I3[) this power determines the 
magnetic field noise by 



6B(f) = VSbJJ) = 



V^kT^PU) 



ANIuj 



(4) 



Since the power P scales quadratically with the driving 
dipole p = AN I in the linear response regime, the above 
equation is independent of the size and driving current 
of the pickup coil. The usefulness of this expression lies 



in the fact that in most cases, calculation of power loss is 
much easier than that of magnetic field noise, the latter 
requiring incoherent sum of vectorial contributions from 
many fluctuation modes inside the source material. 

For high-permeability metals and ceramics used for 
magnetic shields the primary sources of power loss at 
low frequencies (< 1 MHz) are eddy current loss Poddy = 
^ y i^aE 2 dv and hysteresis loss Ph ys t = J v ' H 2 dv 



15j . Here a is the conductivity, p!' is the imaginary 
part of the permeability p — p! — ip," , and the integrals 
are over the volume of the material in which oscillating 
electric and magnetic fields of amplitude E and H, re- 
spectively, are induced by I(t). 

For a given driving dipole strength p, the eddy cur- 
rent j = aE is proportional to the frequency u>, therefore 
Poddy leads to a frequency independent (white) noise ac- 
cording to Eq. to the extent that a is frequency inde- 
pendent. On the other hand, Ph ys t, assuming frequency- 
independent p, leads to a noise with 1// power spectrum, 
which is indeed observed in experiments with ferromag- 
netic transformer cores (T3 |. In what follows we will de- 
note the noises associated with P e ddy and Ph ys t by SB curl 
and SB magn , respectively. 



III. POWER LOSS CALCULATION FOR 
HIGH-PERMEABILITY SHIELDS WITH 
CYLINDRICAL SYMMETRY 

In this section we calculate power dissipation in high- 
permeability shields with cylindrical symmetry when the 
driving dipole is on and along the axis of the shield. See 
Fig. 1(a) for a representative geometry. We restrict our- 
selves to quasi-static regime where the magnetic field 
amplitude inside the shield material is given by its dc 
value, ignoring perturbation due to induced (eddy) cur- 
rents which is proportional to the frequency. The power 
dissipation when the dipole is at other locations and 
along other directions can be calculated numerically with, 
for instance, a three-dimensional finite clcmcmt analysis 
software commonly used for power loss calculations in 
transformer cores. 

Figure 1 (a) also shows several magnetic field lines, cal- 
culated numerically, in the p-z plane around the shield 
generated by a current loop modeling a driving dipole. 
Two features are noticeable. First, the field lines en- 
tering the shield are very nearly normal to the surface, 
reflecting the well-known boundary condition involving 
a high-permeability material. Second, most of the field 
lines reaching the shield arc subsequently confined within 
the thickness of the shell, running nearly parallel to the 
profile of the shield in the p-z plane. This, combined 
with the condition V x B = 0, requires that the field 
lines are nearly uniformly spread within the thickness of 
the shield. For a shield surface with radius of curvature 
(in the p-z plane) R Cl it can be shown that the variation 
of the field strength across the thickness t of the shield is 
6Bh/B» w t/R c , where Bh is the field component parallel 
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FIG. 1: (a) Magnetic field lines inside a high-permeability 
(fir ~ 1000) shield with cylindrical symmetry (mirror sym- 
metry is not assumed). The field is generated by a current 
loop, carrying 1 A dc, whose cross section is indicated by a 
small circle. The eleven field lines enclose magnetic flux of 
nA$,n = 1,2, ••■ ,11, where A$ = 10 -10 Wb. (b) Cross 
section of the shield in the p-z plane, (c) Four geometries 
considered in the power loss calculation: infinite plate, sphere, 
infinite cylinder, and finite closed cylinder. The cylinders and 
the sphere are hollow shells of thickness t 



to the shield in the p-z plane. The condition for the field 
confinement can be estimated, from dimensional consid- 
eration, to be p r t/a 1, where p r is the relative perme- 
ability, and a is the characteristic distance between the 
driving dipole and the shield surface [l6| . Since the same 
factor p r t/a also determines the shielding factor we 
can assume this condition is satisfied if the shell is to 
function as a magnetic shield in the first place. In sum- 
mary, we assume the following for our calculations: (1) 
p r 3> 1 so that the normal entrance boundary condition 
is satisfied. (2) p r t/a ^ 1 so that most of the field lines, 
once entering the shield material, arc confined within the 
thickness of the shield. (3) t/R c <C 1 for most part of the 
shield so that the confined field amplitude is uniform in 
the direction normal to the shield surface [l7j . 



A. Eddy-current loss 

Suppose that the driving dipole is oscillating sinu- 
soidally at a frequency to, p(t) — pzsmujt. We want 
to calculate, to the lowest order in oj, the eddy current 
in the shield which is symmetric around the z-axis. We 
define the position of an arbitrary point in the shield 
in the p-z plane by coordinate (I, s) as shown in Fig. 
1(b). Here I defines a position in the midplane of the 
shield by measuring its distance from the z axis along 



the cross-section of the shield. The coordinate s repre- 
sents the normal distance of a point from the midplane, 
— t/2 < s < t/2. Since we are interested in a thin- walled 
shell, we ignore the variation of the radial coordinate p 
on s: p{l,s) « p(l, s = 0) = p(l). Our assumptions in 
the preceding paragraph imply that the magnetic field 
within the shield material is parallel to the line defining 
the / coordinate, and its amplitude B\\ = B\\ (I) depends 
only on I. Finally we define B±(l) as the amplitude of 
the magnetic field entering the inner surface of the shield 
at (l,s= -t/2). 

In three dimensions, a point (I, s) corresponds to a ring, 
and we define s) as the amplitude of the flux gen- 
erated by the driving dipole p(t) that threads the ring. 
Then the amplitude of the eddy current flowing along the 
ring is 



j = oE^ = cr • lu ■ s)/2np(l). 



(5) 



If all the field lines are confined within the shield, a ring 
on the outside surface of the shield has no net flux in it, 
$(Z, s = t/2) = 0. For all other s, V • B = dictates that 



<!>(/. .si = 2^(7) /)(/). 



From Eqs. (|5l6p the eddy-current loss is 

,! m „ ft/1 ^ 



(6) 



eddy 
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-s) Bi(l)2irp(l)dsdl 



yrn^j- 1 BfAl)p(l)dl 
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1 



6/2 



t 

S 

2 



ds 



^nauj 2 tp, 

where the configuration integral f3, having a dimension 
of flux squared, is 



(3= / t 2 BHl)p(l)dl. 
Jo 



(7) 



This expression can be reduced to a form more useful 
in practical calculations by expressing B\i (I) in terms of 
B ± (l). From V • B = it follows that tB\\(l)p(l) = 
J^B ± (l')p(l')dl'. Therefore, 



= 



B ± {l')p{l')dV 



B. Hysteresis loss 



The hysteresis loss arises from a phase delay in the 
magnetic response of a material to the applied oscillating 
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magnetic field. For most soft magnetic materials used for 
magnetic shields, this delay is small at frequencies below 
~1 MHz. In the following we assume that the shield 
has a constant permeability throughout its volume with 
fi" <C fi' ~ /i r /io- The expression for Ph ys t, to the first 
order in fj,", can then be obtained as follows. 



P, 



hyst 



up," H 2 dv 



l ' 2 1 p" 
-uj — 

ii J-t/2 2 M 



uj!--rB 2 (l)2irp(l)dsdl 



Therefore, both P e ddy and Physt are proportional to (3. 
It follows that the ratio between magnetization- and 
current- induced noises in a cylindrically symmetric shell 
measured on and along the axis is 



5B„ 



Physt\^_ / 3/i" 



1/2 



3 (5 s kin rr j 

V tan di oss , 



SB CUII \P cdd J {au/jPty V2 t 

(9) 

where we used the definitions of skin depth <J s kin = 
1 / y/ir fi' a f and loss tangent tan<5i oss = //'///. There- 
fore 5-Bmagn becomes relatively important when the skin 
depth is greater than ~ t/^tan <5i oss . This is equivalent 
to / < /magn where 



/magn = 3 tan <5l OS s / 2lT p r jlQat 



C. Field noise equations 



(10) 



In this section we list explicit formulas for the magnetic 
field noise for shields of simple geometries shown in Fig. 
1(c), namely an infinite plate, infinite cylindrical shell, 
spherical shell and a finite-length, closed cylindrical shell. 
From the considerations in the previous sections, the on- 
axis mangetic field noise inside a cylindrically symmetric, 
thin- walled shield can be calculated analytically from the 
knowledge of B±(l). Calculation of B± is analogous to 
that of an electric field on the inside surface of a con- 
ducting shell induced by an on-axis electric dipole. Such 
calculation is most easily performed by the method of 
an image in case of an infinite plate and a sphere. For 
a cylinder, Smythe [l8| gives a series expansion solution 
that can be readily adopted for calculation of B±. 

infinite plate 

The midplanc of the plate is the x-y plane, and the 
driving dipole pz is at z = a on the z-axis. I is mea- 
sured from the origin. Due to the image effect of a high- 
permeability plate, B±(l) is twice as large as the normal 
component of a dipolar field expected in free space. Ex- 
plicitly, 



Hop —1 + 3 cos 2 
2?r (a 2 + l 2 f> 2 



where cos 6 



a/y/a 2 
SB cun 



I 2 . This gives 



1 HoVkTat 
/6n a 



(11) 



spherical shell 

For a driving dipole pz at the center of a sphere with 
radius a, the image "dipole" consists of two "monopoles" 
±2ap/d 2 positioned at z — ^f2a 2 /d, in the limit d — > 0. 
The resulting surface normal field is 



B±(l) 



Sfi p I 
- — ^cos- 



where / runs from the north pole to the south pole of the 
sphere, < / < 7ra, and 



SB C 



1 Ho\/kTat 
/2?r a 



(12) 



infinite cylindrical shell 

Smythe [18| gives the electrostatic potential V(p, z) in- 
side an infinitely long conducting cylindrical tube, sym- 
metric around the z-axis, due to a point charge q in- 
side the tube. When q is at the origin and the tube is 
grounded, it is 



V{p,z) = 



(I 



2TT€ a 



a J (ap/a) 
U 2 {a) ■ 



where eo is the permittivity of vacuum, J n (x) is the 
Bessel function of order n, and the summation is over 
the zeros of J ; Jo (a) = 0. From this expression, the 
surface normal (radial) magnetic field at p = a due to a 
magnetic dipole pi at the origin can be obtained as 



B±(z) 
As a result, 

<5P C UIT = 

G = 



d d 

»oead-Q^-Q-V(p,z)\ p = a , qd=v 



sign(z) 



2?ra 3 



*|z|/a 



■A (a) 



/'o 



VkTcrt l~2 



3tt 



G. 



(13) 



Ji(a) 



0.435. 



closed cylindrical shell of finite length 

Ref. [18| also gives the electrostatic potential when the 
conducting cylinder is closed, at, say, z = ±L/2, by con- 
ducting plates . For a charge q at (p = 0, z = Zi), the 
potential at a point (p < a, z > z±) is 



V(p,z; zi) = 

q 



^ sinha(^ 



^l) sinha(^ 

a / V 2a 



J (ap/a) 



sinh a - 



a J 2 (a) 
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If the conducting shell is replaced by a high-permeability 
magnetic shield and a magnetic dipole pi replaces q, the 
normal magnetic field at the top plate is 



B 



top 



B z (p, z = L/2) 
d d 

-fj. eod——V{p,z; zi)\ z = L / 2 , q d= P - 



dz dz\ 

Similarly the normal field on the side wall at z > z\ is 

B s l de = B p (p = a,z) 
d d 

= -Haeod——V(p,z;z 1 )\ p=atqd= p. 

For simplicity, in the following we consider only the case 
when pi is located at the origin, z\ = 0, which gives 
the noise at the center of the shield. Then by symme- 
try calculation of (3 requires integral over only the up- 
per half of the cylinder. The integral path consists of 
two portions: the top plate where / runs along the line 
(0 < p < a, z — L/2) and the upper half of the side wall 
where I runs along the line (p = a, L/2 > z > 0). 
Explicitly, 



1 [ a 1 
^= / -dp 

2 Jo P 



BT(p')p'dp' 



L/2 



1 



-dz 



jy± P (p') P' d p' + J^B^iz') adz' 



The first term can be evaluated using Bessel 
function identities Jq u' Jq (v!) du 1 = u.Ji(u) and 

J Q dxx Ji(ax)Ji(a'x) — jJ 2 (a)5 aa i . This turns out to 
^ (^) 2 2^1 (L/a), where 



1 



sinh 2 # J!(a)' 



(14) 



The second term is more tedious, but can be reduced 
to (^) 2 ^F 2 (L/a) with 

E 



1/2 L 

dx— 



cosh 



aLx 



Ji(a) 



sinhf^ 

2 a 



(15) 



Finally the field noise is 



5B r 



pLo^JkTat I 2 



3?r 



G, 



(16) 



G = F 1 (L/a) + 2F 2 (L/a). 



Numerical evaluation of the above equation shows that 
G = 0.657, 0.460, 0.438 for aspect ratios L/2a = 
1, 1.5, 2, respectively. Thus the noise from a closed cylin- 
drical shield with aspect ratio of 2 already approaches 
that of an infinitely long shield within 0.5%. 



IV. COMPARISON WITH NOISE FROM 
NONMAGNETIC CONDUCTING SHELLS 

An interesting question is how the magnetic field noise 
in a high-permeability shield compares with that in a 
non-magnetic shell with the same geometry and con- 
ductivity. As indicated in Ref. [4|, calculation of low- 
frequency eddy current loss in an axially symmetric, non- 
magnetic metal driven by an axial dipole p — pi sin tot is 
relatively simple, because the amplitude of the induced 
electric field is proportional to the magnetostatic vector 
potential (in Coulomb gauge) due to a dipole in vac- 
uum. For an axial dipole pi at the origin 



A<j,{p,z) = 



PoP 



4tt 



(p 2 + z 2 ) 



3/2 



and 



-Poddy = 



Aidv. 



Equations for the quasi-static field noises associated 
with this loss are listed in Table |T] for the geometries 
considered in the previous section. It is found that the 
current-induced noise inside a high-permeability shell is 
in general not much different from that inside a non- 
magnetic shell. The difference can be either positive (in- 
finite plate) or negative (sphere and cylinder). Qualita- 
tively, one can think of two competing effects, namely 
self-shielding and image effects, due to the high perme- 
ability of the material. In a long tube, the field generated 
by a noise current at the end of the tube is self-shielded as 
it propagates inward. On the other hand, the field gener- 
ated by a current loop on the surface of an infinite plate 
is amplified because of an image current adding field in 
the same direction. 

A dramatic illustration of the latter effect is found 
in the case of field noise in between two infinite plates, 
with thickness t, separated by L. When the plates are 
non-magnetic, the total quasi-static power loss induced 
by an axial driving dipole half way between the plates 
is simply twice that induced in a single plate. In the 
limit pirt/L — > co, however, it can be shown that the 
power loss and therefore the noise logarithmically di- 
verges. This is because the noise current in either plate 
generates an infinite series of image currents, and when 
all the current modes are considered their contributions 
do not converge. It is evident that the noise from a high- 
permeability structure, even in the quasi-static regime, 
cannot in general be obtained from the quadrature sum 
of the noise from its individual parts. 



V. FREQUENCY DEPENDENCE 

Here we consider how the noise 5B cmv considered in 
Section III and IV rolls off at frequencies above the 
quasi-static regime. Previous theoretical and experimen- 
tal works on noise from conducting plates and enclosures 
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geometry 


field noise due to Johnson noise current 


high-permeability 


non-magnetic 


infinite plate 


1 ^oVkTat 

£>±>curr = , 

V67r a 


_ 1 fioVkTat 
V8tt a 


spherical shell 


1 ^oVkTat 

O-Dcurr = ; 

\/2tv a 


SB _ 12 voVkTat 
V 37r a 


infinite cylindrical shell 


jcd [2G fioVkTat 

0£>curr=\ /- , G ^0.435 

V 6tv a 


_ /~3~ fiaVkTat 
V 16 a 


finite, closed 
cylindrical shell 


Eas. (|14I15I16P 


6B = VG IM,y/hT,Tt , 
a 

1 /3(L/2a) 5 +5(L/2a) 3 + 2 ,3,-!^ 
8tt V (i/2a) 2 (1 + (L/2a) 2 ) 2 2a ^ 



TABLE I: Magnetic field noise from high-permeability and non-magnetic plate and shells of conductivity a. The geometries 
are shown in Fig. 1(c). 



references 


frequency dependence 


material 


method 


Rcf. [11], Eq. (5) 


f - /- 1 - r 3/4 


non- or weakly magnetic slab 


calculated 


Ref. [3], Fig. (6) 


f° - r 1 - r 3/4 


nonmagnetic slab 


calculated 


Ref. [7], Fig. (1) 




high-permeability slab 


calculated 


Ref. [9], Fig. (2) 




nonmagnetic, thin sheet 


calculated 


Ref. [20], Fig. (2) 


/° - r 1/4 


mu-metal plate 


measured 


/° - r 1 


copper plate 


measured 



TABLE II: Frequency dependence of magnetic field noise induced by Johnson noise current in magnetic and non-magnetic 
metallic plates. 

[20| reported initial roll-off given by SB(f) oc f~ J , where 
7 w 1 for non-magnetic metals and 7 w 1/4 for high- 
permeability metals. Below we provide qualitative ex- 
planation of such dependences by considering a simple 
model. 

Suppose we measure noise from a large, thin plate with 
conductivity a at a distance a along the direction perpen- 
dicular to the plate. We assume that a is independent of 
frequency. The plate has a thickness t <§C a and a lateral 
dimension much larger than a. It is reasonable to as- 
sume that the field noise mostly comes from fluctuating 
currents flowing in a series of concentric rings directly 
below the measurement point with radius on the order 
of a. Since these current paths are connected in parallel, 
we can assume that in fact the noise comes from cur- 
rent fluctuation in a single annular loop of mean radius 
« a and width w a. The dc resistance of such a loop 
is i?o = 2tt/ at, which gives conventional Johnson noise 
current SI = ^AkT/R Q « y/{2/n)kfat. The magnetic 
field noise arising from this current is indeed of the same 
order of magnitude as the noise calculated in the previous 
sections. 



the skin depth effect to R(f > / s kin) ~ 27r/cr^ s km oc f 1 ^ 2 . 
The threshold frequency is 

/skin = l/7TMrM0Ci 2 . (17) 

Second, the self inductance L of the loop suppresses SI 
if 2irfL > R(f). Therefore the current noise should in 
general be written as open-loop voltage noise divided by 
total impedance, 

= ^4kTR(f) 
\R(f) + i2nfL\ ■ 

where R(f) includes the skin depth effect. If the condi- 
tion 2ttJL > R(f) is reached at a frequency f- ln d < /skin, 
such frequency is obtained from 2irf m dL = Rq, namely, 

f ind = 1/atL = 1/Cfioata, (18) 

where C is a constant of order unity. For a non-magnetic 
plate, /ind//skin=(7r/C)/x r t/a <c 1 and inductive screen- 
ing indeed appears at a frequency far below that at which 
skin depth becomes important. The initial roll-off of the 
noise then occurs at / > /md, where the current noise 
scales with frequency as 



At high frequencies this current is suppressed in two 
ways. First, when <5 s ki n < t, the resistance increases by 



7 



As / further increases beyond / s km, the scaling changes 
to 



/skin ^ f- 



On the other hand, for a high-permeability plate 
used for magnetic shields, skin depth effect appears 
at a frequency far below that for inductive screening, 
/ind//skin 3> 1. Therefore the initial roll-off is expected 
to follow 



SI: 



R(f) 



-1/4 



/skin < / < //, 



ind ■ 



The frequency f( nd at which inductive screening becomes 
important for a high-permeability plate is obtained from 
2n fLd L = R (fLd) = 27r/cr(5 s kin which reduces to 



fLd = (tt/C 2 ) Hr/^ocra 2 . 



(19) 



Beyond this frequency (5/ again scales as /~ 3 / 4 . 

Table [TT] summarizes the frequency dependence of the 
Johnson-current-induced magnetic field noise reported 
in five references. It is found that our simple model 
correctly predicts all the essential features of the fre- 
quency dependences found in these works. For non- 
magnetic plates, the two threshold frequencies Eq. (fTT]) 
and Eq. (JTSJ) agree, up to a numerical factor, with those 
obtained in Ref. [| [U and Ref. [01 H- For high- 
permeability plates, Table (1) of Ref. Q also can be inter- 
preted as giving the same threshold frequencies between 
different regimes, Eq. (fT7|) and Eq. (fT9| . obtained in this 
work HI]. 

Finally, if we include the magnetization-fluctuation 
noise calculated in Section III-B, the magnetic field noise 
from a high-permeability plate is expected to exhibit a 
rather complicated frequency dependence 



SB(f) : / 



-1/2 



-1/4 



/ 



-3/4 



where the three threshold frequencies dividing different 
scaling regimes are given by Eq. (fTOjl . Eq. (fT7)) , and 
Eq. (|19|) . in the increasing order. 



VI. NOISE REDUCTION BY DIFFERENTIAL 
MEASUREMENT 

A common technique to reduce the effect of magnetic 
field noise from a distant source is to make a differential 
or gradiometric measurement. In the first-order differen- 
tial measurement, one measures Sdiff(i) — B\{t) — B%(t), 
where B\ and Bi are the magnetic fields at two points 
separated by a baseline d. The fluctuation in this quan- 
tity SBdis(f) can be calculated following the same prin- 
ciples described in Section II, with a single pickup coil re- 
placed by two coils connected in series so that the induced 
voltage is proportional to -Bdiff- Reduction of noise from 
a distant source now corresponds to reduction of power 



loss induced in the material when driven by this "gra- 
diometric" coil, which appears as a quadrupole, rather 
than a dipole, seen from a distance a ^S> d. If the two 
coils connected in series are identical, each represented 
by an oscillating dipole of amplitude p, than the result- 
ing power loss P gives <5-Bdiff through 



SB dm (f) 



V4fcTv/2P(7) 



In the limit a > rf, P is proportional to the square of 
the driving quadrupole moment p 2 d 2 . From dimensional 
consideration, therefore, 5B<uh scales as (d/a). 

Table [nil shows the results of analytical calculations of 
SBdift for an infinite plate and an infinitely long cylindri- 
cal shell. Only the white noise associated with the eddy 
current loss is considered. The noise is calculated for an 
axial differential measurement along the symmetry axis, 
in the limit where the baseline is much smaller than the 
shortest distance a to the material. It is seen that in all 
cases the noise reduction factor is very nearly d/a. 



VII. CONCLUSION 

We have used generalized Nyquist relation applied 
to electromagnetic power dissipation and magnetic field 
fluctuation to calculate magnetic field noise inside high- 
permeability magnetic shields. Analytical results for axi- 
ally symmetric geometries show that the quasi-static field 
noise due to Johnson current noise in a metallic shell is 
slightly altered as the material gains high magnetic per- 
meability. For magnetic shields with small electrical con- 
ductivity, 1// noise from magnetization fluctuations be- 
comes dominant over Johnson-current-induced noise be- 
low a threshold frequency proportional to its magnetic 
loss factor. Established numerical methods of finite- 
element analysis of electromagnetic power loss can be of 
great utility in calculating magnetic field noise spectrum 
from dissipative materials of complicated geometry. At 
relatively high frequencies, one could experimentally de- 
termine the power loss in dissipative materials using a 
pickup coil. This has an advantage that no prior knowl- 
edge of material parameters is necessary to predict the 
field noise. From Eq. ([3]) and ([4|, it turns out that a 1 
fT/Hz 1 / 2 noise at 1 kHz and at room temperature corre- 
sponds to an effective resistance of 10 mil in a 1000-turn 
driving coil of 5 cm diameter. This change in the re- 
sistive load is within the measurement range of modern 
impedance analyzers. 

As reported earlier [2J, we find that quasi-static John- 
son current noise in magnetic shields is significantly 
higher than intrinsic noise of modern magnetometers. 
Due to a small skin depth of high-permeability materials, 
however, the white noise range extends only to relatively 
low frequencies (/skin = 1~T00 Hz), beyond which the 
noise rolls off as / _1 / 4 , until self-induction effect fur- 
ther brings down the noise. This indicates that usual 
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geometry 


material 


SB diB 




infinite plate 


high-permeability 


1 HovkTat d 
Vin a a, 


1.22* 
a 






3 


non-magnetic 


/ 3 fioVkTat d* 
V 167T a a 


1.22* 

a 


s \ 






t 


infinite cylindrical shell 


high-permeability 


j~2G fj, y/kTcrt d 
V 37r a a' 

J.J Z Ji(a) ) abi8 


1.19 d - 

a 




A 


i 


non-magnetic 


/ 45 fj-oVkTat d 
V 256 a a 


0.97^ 

a 







*This agrees with Eq. (43) of Ref. | 

TABLE III: Differential measurement noise from an infinite plate and a long cylindrical shell of conductivity a. Arrows on two 
circular loops in each of the inset figures indicate the direction of the magnetic fields contributing to the measurement. <5-B S m g ie 
in the last column is the magnetic field noise in non-differential measurement taken from Table [I] 



geometry 


SB CUII [fT/Hz 1 / 2 ] 


5B m&m [fT/Hz 1 / 2 ] at 1 Hz 


numerical 


analytical 


numerical 


analytical 


infinite plate 


3.63 


3.68 


2.01 


2.07 


spherical shell 


6.38 


6.38 


3.57 


3.59 


closed cylindrical shell 










L/2a = 0.5 


12.4 


12.2 


6.93 


6.87 


L/2a = 1.0 


6.01 


5.97 


3.33 


3.36 


L/2a = 1.5 


5.04 


4.99 


2.77 


2.81 


L/2a = 2.0 


4.92 


4.87 


2.70 


2.74 



TABLE IV: Magnetic field noise calculated for mu-metal plate and enclosures with a = 1.6 x 10 6 f2 _1 m _1 , fi r — 30,000, 
tan 5 = 0.04. Geometrical parameters are a = 0.2 m, t = 1 mm, referenced to Fig. 1(c). Column 3 is calculated from equations 
in column 2 of Table 1. Column 5 equals column 3 multiplied by 0.5628, from Eq. ([9}. Numerical calculation for an infinite 
plate was obtained by extrapolation of the results for finite-size plates. 



geometry 


^^[fT/Hz 1 / 2 ] 


numerical 


analytical 


infinite plate 


15.4 


15.5 


spherical shell 


35.7 


35.9 


closed cylindrical shell 






L/2a = 0.5 


45.0 


44.9 


L/2a = 1.0 


34.1 


34.2 


L/2a = 1.5 


33.6 


33.7 


L/2a = 2.0 


33.6 


33.7 



TABLE V: Magnetic field noise from eddy current loss cal- 



culated in aluminum, <7 = 3.8xl0 7 f2 
geometries are the same as in Table HVl 



fj r 



1. The 



room-temperature mu-metal shields may be used with- 
out adding significant noise if the signal is modulated 
at relatively high frequencies. At low frequencies, most 
sensitive experiments would require a low-loss noncon- 



ducting magnetic materials, such as certain ferrites, as 
the innermost layer of a multi-layer shield, or differential 
field measurement with a short baseline. In practice, a 
combination of these techniques should be implemented 
to suppress shield-contributed noise to an insignificant 
level. 



APPENDIX A: COMPARISON WITH 
NUMERICAL CALCULATIONS 

Here we compare magnetic field noises predicted by an- 
alytical expressions in Table U with those obtained from 
numerical calculations of power loss for representative 
geometries. The calculation was performed by a finite 
element analysis software (Maxwell 2D, Ansoft) which 
determined electromagnetic fields in space on a mesh 
through iterative solution of the Maxwell's equations. 
The driving dipole was modeled as a small current loop 
on the symmetry axis. For P e ddyi the current oscillated 
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(a) 



(b) 



(c) 



& 



(d) 
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FIG. 2: Definition of geometries used to calculate magnetic 
field noises from (a) a small solid sphere, (b) a thin disk with 
arbitrary diameter, (c) a planar, infinite thin film divided into 
a square array of small disks, and (d) a long, thin wire with a 
circular cross-section. In (c), the small disks are close-packed, 
but electrically isolated from each other. 



at / = 0.01 Hz. For Physt, a magnetostatic problem was 
solved with a static current in the coil, and volume inte- 
gral of H 2 in the material was calculated. Magnetic field 
noises were then obtained by Eq. ([4]) . The errors due to a 
non-zero radius of the loop were insignificant within the 
accuracy of the numerical calculations presented here. 

Table IIVI shows magnetic field noises from high- 
permeability plate and shields. The loss tangent assumed 
is for illustration purpose only. It is seen that in all cases 
considered here, numerical and analytical results differ 
by less than 3%. The errors represent the accuracy of 
the assumptions made in magnetic field calculations in 
Section III. 

Table[V] s hows magnetic field noises from non-magnetic 
plate and shells. These numbers can be used to estimate 
noises from non-magnetic, metallic enclosures often used 
for radio-frequency shielding. The differences between 
analytical and numerical calculations, less than 1%, are 
consistent with the errors in the numerical calculations. 



was calculated by the method outlined in section IV. For 
others, the eddy current density can be calculated from 
the equations 



V x j = crV x E 
V ■ j = 



with the boundary condition that the normal component 
of j is zero on the surface (boundary) of the object. Here 
B is the amplitude of the oscillating magnetic field gen- 
erated by p(t) in free space. For a thin film lying in the 
x-y plane, V x j is along the z axis and therefore only 
B z contributes to the loss. The two-dimensional current 
distribution j(x, y) = (u(x,y),v(x,y)) then satisfies 



du dv 

dy dx 

du dv 

dx dy 



iaujB z 
0. 



(Bl) 
(B2) 



When the film is divided into small patches whose lat- 
eral dimensions are much smaller than the distance to 
the dipole, the current distribution in each patch can be 
calculated assuming a constant B z within the patch. 

When the film is in the shape of a long, narrow strip, 
such as a long straight wire patterned on an insulating 
substrate, the noise measured along the z axis on a point 
in the x-y plane can be calculated by solving Eqs. (|B1[ 
IB2[) with the boundary condition u(x = ±L/2, y) = 
0, v{x, y = ±j/ ) = 0. Here the strip is assumed to 
occupy a region —L/2 < x < L/2, — yo < y < yo 
with L 3> yo- The source term is given by B z (x,y) = 
(/j,op/47ra 3 )(l + x 2 /a 2 )~ 3 / 2 , assuming the noise is mea- 
sured at (0, a, 0) and a 3> yo- Eqs. (|B11 IB2[) are then 
satisfied by 



u(x,y) 
v(x,y) 



aujfiop ^ sinh k n y 



iraL 

n 

naL 2— 1 



sin k n x — 1 



cosh k n yo 

cosh k n y 



cosh k n y Q 



Ki{ak n ) 



APPENDIX B: MAGNETIC FIELD NOISE FROM 
OTHER METALLIC OBJECTS 



where K\ is the modified Bessel function of order one 
and k n = (2n — l)ir/L, n = 1,2, • • • . The power loss in a 
strip with thickness dt, calculated in the limit L/a — > oo, 



For the purpose of future reference, here we list equa- 
tions for magnetic field noises resulting from Johnson 
noise currents in non-magnetic, conducting objects with 
simple geometry. We only consider white noise in the 
low frequency limit. Table IVT1 lists equations for a small 
solid sphere, thin planar films, and a long thin wire, as 
defined in Fig. @ . In the context of an atomic vapor-cell 
magnetometer, these objects can be associated with an 
alkali metal droplet, low-emissivity conductive coatings 
on a glass, and a heating wire, respectively. For prob- 
lems with a cylindrical symmetry (Fig. ((2^,b) ) , the eddy- 
current loss induced by a driving dipole p(t) = psinwi 



P(y )dt = - 
a 



(u 2 + v 2 ) dxdy 



For a wire with a circular cross-section, the integral of 
the above equation over the profile (in the y-z plane) of 
the wire gives the total loss 



p . — 



f_° yo p (^y*- t2 ) dt - 
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geometry 


SB 


figure 


small solid sphere 




(2a) 


thin film disk 


i f-io \ Ki at i 


(2b) 


^ a 1 + a 2 /r 2 


infinite square array of small disks 


1 ^ n n \/hTfTt 7 
V 2048 a 2 


(2c) 


circular cross-section wire 




(2d) 



TABLE VI: Magnetic field noise due to Johnson noise current in metallic objects with conductivity a. 



The corresponding magnetic field noise, in the leading 
order in yo/a, is 
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SB 



3 

128 



5/2 



For example, a long straight constantan wire with di- 
ameter 2yo = 1 mm and a = 2xl0 6 f2 -1 m -1 exhibits 
5B = 0.433 fT/Hz 1 / 2 at room temperature when mea- 
sured at a = 1 cm in the direction perpendicular to both 
the wire and the normal direction to the wire. 
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